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ABSTRACT 

Extra-galactic radio sources are a significant contaminant in cosmic microwave 
background and Sunyaev-Zel'dovich effect experiments. Deep interferometric ob- 
servations with the BIMA and OVRO arrays are used to characterize the spatial, 
spectral, and flux distributions of radio sources toward massive galaxy clusters 
at 28.5 GHz. We compute counts of mJy source fluxes from 89 fields centered 
on known massive galaxy clusters and 8 non-cluster fields. We find that source 
counts in the inner regions of the cluster fields (within 0.5 arcmin of the cluster 
center) are a factor of 8.91 4 /! times higher than counts in the outer regions of 
the cluster fields (radius greater than 0.5 arcmin). Counts in the outer regions 
of the cluster fields are in turn a factor of 3.3+ 4 § greater than those in the non- 
cluster fields. Counts in the non-cluster fields are consistent with extrapolations 
from the results of other surveys. We compute spectral indices of mJy sources 
in cluster fields between 1.4 and 28.5 GHz and find a mean spectral index of 
a = 0.66 with an rms dispersion of 0.36, where flux S oc v~ a . The distribution 
is skewed, with a median spectral index of 0.72 and 25th and 75th percentiles of 
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0.51 and 0.92, respectively. This is steeper than the spectral indices of stronger 
field sources measured by other surveys. 

Subject headings: galaxies: clusters — cosmology: observations — cosmic mi- 
crowave background — radio continuum: galaxies 



INTRODUCTION 



Extra-galactic radio sources are a significant contaminant in cosmic microwav e back - 
ground (CMB) and Sunyaev-Zerdovich effect (SZE) experiments (e.g., see iHqlder J2002h: 

Knox et_ali J2004h : iTegmark et all (boooh ) . Measurements of the CMB and SZE JSunvaev k Zel'dovich 
1970l . Il972h have the potential to yield a wealth of cosmological information if foreground 
contaminants are well-understood. Models for the number counts of radio sources as a 



function of flux (IDe Zotti et al.ll2005l ; iToffolatti et al.lll999t ISokasian et al.ll200ll ) have been 
derived from observations at lower frequencies and extrapolated to microwave frequencies. 

Radio sources are often associated with the clusters of galaxies themselves. This is a 
potential source of bias for current and planned SZE surveys such as the Sunyaev-Zel'dovich 
Array the Arcminute Microkelvin Imager (AMljj, the Atacama Pathfinder Exper- 

iment Sunyaev-Zel'dovich (APEX-SZ) survejo, the South Pole Telescope and the 

Atacama Cosmology Telescope (ACT)|. To understand the impact on planned SZE cluster 
surveys, it is critical to characterize the spatial, spectral, and flux distribution of sources 
associated with clusters. 

Contaminating radio emission from extra-galactic sources at frequencies less than ap- 
proximately 100 GHz is attributed to synchrotron radiation from active galactic nuclei (AGN) 
and star-forming galaxies. The AGN-powered radio galaxies dominate the source counts at 
high luminosities. At higher frequencies, the emission is attributed to dust emission from 
star-forming galaxies. ( 



Slee et al. 


(1983. 


1998 


); 


Owen 


(1996 


); 



Reddy fc Yunl (120041 )) show a strong central concentration of radio galaxies in clusters. The 



*SZA website: http://astro.uchicago.edu/sza 
2 AMI website: http://www.mrao.cam.ac.uk/telescopes/ami/ 
3 APEX-SZ website: |http : / /bolo .berkeley. edu / apexsz 



4 SPT website: |http: / /spt.uchicago.edu/ 



ACT website: http://www.hep.upenn.edu/act/ 
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distribution of synchrotron-em itting; star-forming galaxies is found to be less centrally-peaked 
in clusters (IRizza et al.ll2003l ) than that of AGN-powered radio galaxies. At our sensitivity 
level and observing frequency of 28.5 GHz, radio sources powered by AGN dominate and we 
refer to them henceforth as radio sources. 



Several CMB experiments such as W MAP (Bennett et all 120031). DASI flKovac et al. 



20021 ) . VSA (jCleary et al.ll2005l ). and CBI (IMason et al.l 120031 ). have measured microwave 
source counts as a function of flux for sources brighter than about 10 mJy. There have 
been two prior analyses of radio sources in SZE data taken with the Owens Valley Radio 
Observatory (OVRofl and Berkeley-Illinois-Maryland Associ ation (BIMAfl array s. Using a 
sample of 56 fields centered on known massive galaxy clusters, ICooray et al.l (119981 ) computed 
counts and spe c tral i ndices of radio sources. Using the outer regions of 41 cluster fields, 
LaRoque et al.l (120021 ) computed the normaliza tion of source count s as a function of flux 



for fai nt sources in SZE data. The data used in ICooray et al.l (119981 ) and in lLaRoque et al. 
( 120021 ) are subsets of the data presented in this paper. 



In this paper we analyse faint (~ mJy) radio sources found serendipitously toward 
massive galaxy clusters at 28.5 GHz from the OVRO/BIMA SZE imaging project. Charac- 
terizing the spatial, spectral, and flux distribution at this relatively high frequency should 
help improve projections for radio source contamination in SZE and CMB experiments at 
frequencies of 30 GHz and higher. We use 89 fields centered on known massive galaxy 
clusters and 8 non-cluster fields. The paper is organized as follows: Section 2 reviews the 
observations, data reduction, field selection, and the measured source fluxes. In Section 3, we 
compute spectral indices between 1.4 and 28.5 GHz using fluxes from our data and published 
1.4 GHz surveys. In Section 4, we present source counts as a function of flux for cluster and 
non-cluster fields as well as the angular radial dependence of counts in cluster fields. We 
compare our results with those from other experiments and with theoretical models. 



2. OBSERVATIONS 



2.1. Observations 



The 28.5 GHz observations were carried out with BIMA during the summers of 1996 
2002 and OVRO during the su mmers of 1995 - 20 01 as part of the OVRO /BIMA SZE 



imaging project (see, for example, iReese et al.l (120021 ) and iGrego et al.l (120011 )). A total of 



6 The OVRO mm-wave array is operated by Caltech with support from the National Science Foundation 
7 The BIMA array is operated with support from the National Science Foundation 
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62 cluster fields were observed at BIMA and 55 cluster fields were observed at OVRO. A 
total of 28 cluster fields were observed at both BIMA and OVRO, yieldi ng observations of 89 



uniqu e cluster fields. For the BIMA CMB finescale anisotropy project (IDawson et al.l 12002 



20061 ). a total of 18 non-cluster fields were observed. 



The BIMA array consists of ten 6.1 meter diameter telescopes with primary beams of 
6.6 arcmin FWHM; nine of the ten BIMA telescopes were used for the 28.5 GHz obser- 
vations. The OVRO array consists of six 10.4-meter telescopes, with primary beams of 4.2 
arcmin FWHM; all six OVRO telescopes were used. The primary beams were measured holo- 
graphically and were found to be well-approximated by Gaussians; we use the azimuthally 
averaged measured beam responses in our analyses. For CMB observations at BIMA, the 
array was set in a compact configuration to maximize brightness sensitivity. For cluster 
observations at both OVRO and BIMA, most of the telescopes were also configured in a 
compact configuration that provided dense u — v coverage to the shadowing limit, and one 
or two telescopes were placed at longer baselines for higher angular resolution monitoring of 
sources. The longest baselines used at OVRO and BIMA ranged from 70 to 140 meters for 
the data presented here. 



The telescopes were outfitted with cm- wave receivers ( Carlstrom et al 



9961) equipped 



with cryogenically cooled 26 - 36 GHz HEMT amplifiers (iPospieszalski et al.lll995l ). Typical 
receiver temperatures were 11 K to 20 K, and when integrated with the OVRO and BIMA 
telescopes yielded typical system temperatures scaled to above the atmosphere of 45 K to 
55 K, and as low as 35 K. OVRO observations were correlated with an analog correlator 
consisting of two 1 GHz bandwidth channels centered at 28.5 and 30 GHz. The u — v data 
from the two channels were not combined before making maps or fitting sources to the data. 
The OVRO correlator integration time was four minutes or less. BIMA observations were 
correlated with a multi-channel hybrid correlator. After calibration, the u — v data were 
reduced to a single 0.8 GHz wide bandwidth centered at 28.5 GHz. The BIMA correlator 
integration time was 50 seconds. 

Observations of cluster and CMB fields were interleaved every ~ 20 minutes with obser- 
vations of a strong point source for phase calibration. Source data that were not bracketed 
in time by phase calibrator data were discarded, as were data from baselines in which one of 
the array elements was shadowed by another. Source data were also discarded if the brack- 
eting phase calibrator observations showed a discontinuity in the instrument phase response. 
Lastly, source data were discarded if atmospheric phase fluctuations showed a loss of corre- 
lation greater than a few percent on the long baseline observations of the phase calibrator. 
Observations of the phase calibrators indicate the gain stability was stable to ~ 1% over 
several months. The absolute calibration is based on observations of Mars, with the bright- 
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ness temperature taken from the iRudvl (119871) model. Fur ther details of the observations 
and data reduction can be found in iGrego et al.l (120011 ) and lReese et al.l (120021 ) . 



2.2. Field Selection 

The cluster fields of the OVRO/BIMA SZE imaging project were chosen to obtain 
precise measurements of the SZE in massive galaxy clusters. Potential targets were screened 



for st rong radio sou rces using archival data at lower frequencies such as NVSS ( ICondon et al. 



1998|) and FIRST (IWhite et al.lll997l ). In addition, if a strong source (> 10 - 20 mJy) was 
detected near the cluster center in the initial 28.5 GHz observations, observations ceased in 
favor of other less contaminated targets. Due to the constraints of this field selection, we do 
not attempt to characterize the distribution of sources brighter than 10 mJy in this analysis. 

The cluster fields were chosen mainly from X-ray catalogs and one optical survey, includ 



ing th e following: (1) the RO SAT Brightest Cluster Survey, BCS (lEbeling et al.lll997l . Il998 



2000al ; I Crawford et al.lll999f). (2) the Einstein Observatory Extended Medium Sensitivity 



1994 



Surve y, EMSS (IGioia et al.lll990l : IStocke et al.lll99ll ; lGioia fc Luppino 
1994h . (3) the ROSAT X-ray Brightest Abell Cluste rs. XBACS (libeling et al 



Maccacaro et al. 



(4) the Wide Ang l e ROSAT Pointed Su rvey WARPS (jScharf et al.lll997 



1996b 



Jones et al 



a 



1998 



Fairlev et alil200ol : lEbeling et alJioOObh . (5) the Massive Cluster Survey, MACS (lEbeling et al" 
2001al ). and (6) the Red-Sequence Cluster Survey, RCS (IGladders fc Yed 120051 ) . 



Clusters were selected from the BCS, EMSS, and XBACS X-ray catalogs on the basis of 
high X-ray luminosity. Clusters from WARPS and MACS were selected for X-ray luminosity 
and high redshift. The 8 clusters selected for SZE observations from MACS form a complete 
redshift-selecte d flux-limited X-ray sa mple and were chosen regardless of possible radio source 
contamination (ILaRoque et al.l 120031 ) . A few optically selected clusters were also provided 
by the RCS team. Table 1 lists the coordinates of the pointing centers for the cluster fields. 

The BIMA array was also used to o bserve 18 non-cluster fields for CMB anisotropy 
measurements (jPawson et al.l 120021 . 120061 ) . The pointing center coordi nates for these fields 
are provided in Table 2. Of these fields, only the 8 fields BDF14-BDF21 (toawson et al-lbood ) 
were chosen wi thout regard to pos sible radio source contamination. The fields BDF4, and 
BDF6-BDF13 (jPawson et al.ll2002l ) were chosen from NVSS to have minimal contamination 
from strong radio sources. The field HDF is centered on the Hubble Deep Field and was also 
selected to have no strong radio sources. Only the 8 fields selected without regard to source 
contamination are used in the analysis presented in this paper. 
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2.3. Source Fluxes 

The positions and fluxes of 28.5 GHz sources an d SZE decrements are determined using 



the DIFMAP software package (IPearson et al.lll994j ). The SZE data consist of positions in 



the Fourier domain (also called the u — v plane) and the visibilities — the complex Fourier 
component pairs as functions of u and v , which are the Fourier conjugate variables to right 
ascension and declination. DIFMAP is used to create a map from the u — v data using all 
baselines and natural weighting (oc a~ 2 ). Source and decrement positions are determined 
from this map. Emission sources are initally assumed to be pointlike, while the SZE decre- 
ment (if any) is approximated with an elliptical, isothermal f3 model. The fourier transforms 
of the model components are jointly fit to the observed visibilities. The image of the residuals 
is searched for additional sources at greater than five times the image rms and any sources 
found in this way are added to the model and jointly refitted. Noise estimates are extracted 
from the final residual images. 

The source fluxes determined in this way depend on the coverage of the u — v plane for 
sources with structure on scales larger than that probed by the longest baselines (30 arcsec 
to 15 arcsec for the 70 to 140 meter maximum baselines employed here). To enable a uniform 
comparison with lower frequency data from the NVSS, which has an angular resolution of 
45 arcsec, we repeat the u — v fitting procedure considering only baselines shorter than 4 kX. 
This provides matched angular resolution at a cost of typically 40 — 50% of our data, yielding 
a sensitivity loss of 30 — 40%. We also measured the flux within an extended area centered 
on each radio source in the CLEANed maps made from the 4 kX cut and full u — v data set. 
We compare the fluxes obtained from the u — v model fits and the analysis of the CLEANed 
maps from both the cut and the full u — v data. We find the fluxes found with the four 
methods agree for 103 out of 121 sources, double-counting sources detected with both arrays 
because of the difference in sensitivity and u — v coverage. In these cases we conclude that 
the point source approximation is adequate for our angular resolution and sensitivity and 
we report the fluxes and noise levels determined from all baselines in Tables 1 and 2. The 
remaining 18 sources have best-fit point source fluxes that depend on our choice of u — v 
range. We find, however, that the total flux recovered by the CLEAN algorithm in the region 
around each source in the full data set matches that obtained for the flux found with only 
the short baselines for all but one of these sources. In these cases we use the flux measured in 
the CLEANed map made from the full data set. Finally, in one case (MACS J0717.5+3745, 
source 1) we do not recover the full short-baseline flux when CLEANing the image of the 
full data set. In this case we have used the flux and noise measured for just the baselines 
shorter than 4 kX for both sources in the field. 



Table 1 lists the 28.5 GHz beam attenuation-corrected fluxes for sources detected in 
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the BIMA and OVRO cluster fields and the centroid of the SZE decrement, if any. Noise 
is reported for the field center in column 7, while the uncertainty in the flux of each source 
(column 15) properly accounts for the beam profile. The radial distance from the pointing 
center, or the SZ decrement if observed, to each source is also provided. At the 5<r-level or 
greater we detect 62 sources in the 62 BIMA cluster fields and 55 sources in the 55 OVRO 
cluster fields. A total of 22 sources are detected at both BIMA and OVRO, yielding 95 
unique sources in the 89 cluster fields. For sources observed at both BIMA and OVRO, 
the flux measurements are in good agreement. Source fluxes, positions, and noise levels for 
the non-cluster fields are given in Table 2. We detect two sources at > 5er in the 8 BIMA 
non-cluster fields which were selected without regard to possible source contamination. 

The fluxes in Tables 1 and 2 do not account for attenuation due to temporal and spectral 
averaging of the u — v data, effects which are far less significant than the beam attenuation. 
Averaging of interferometric data in the u — v plane leads to attenuation of the amplitude 
response as a function of angular distance from the field center that depends on the telescope 
array configuration. However, these effects have been restricted to insignificant levels as part 
of the experimental design. The bandwidth constraint translates to ^ x ^ < 1, where ^ 
is the fractional bandwidth, D is the diameter of an array element and B is the baseline. 
For the longest baselines we find values of — x to be 0.32 - 0.64 and 0.24 - 0.47 for the 
BIMA and OVRO observations, respectively. For measurements using only short baselines 
(< AkX) these values decrease to 0.18 and 0.13, respectively. The correspondence between our 
measured fluxes using all baselines and just the short baselines also suggests that bandwidth 
smearing is not important. The integration time constraint translates to ^t int x |j < 1, 
where ti nt is the u — v averaging time in hours. For the longest baselines and integration 
times we find values of ^t int x § to be less than 0.0035 and 0.001 for the BIMA and OVRO 
observations, respectively. 



3. SPECTRAL INDICES 

We use the results of surveys at lower observing frequency to constrain the spectral 
indices of radio sources detected with the BIMA/OVRO observations. Fluxes at 1.4 GHz 
are taken primarily from the NVSS catalog, which has a resolution of 45 arcsec and limiting 
peak source brightness of 2.5 mJy We obtain 1.4 GHz fluxes from the FIRST catalog 
(limiting flux of 1 mJy and resolution of 5 arcsec) for several additional sources which were 
below the NVSS detection threshold. We obtain 1.4 GHz fluxes from VLA archival maps 
for several sources which were not in the NVSS or FIRST catalogs. Of the 95 sources in 
the cluster fields, 88 have unambiguous counterparts in the NVSS, FIRST, or VLA archival 
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data. All 28.5 GHz sources within 0.5 arcmin of the cluster centers were confirmed to have 
counterparts at 1.4 GHz. The 1.4 GHz fluxes and the results of a literature search for low 
frequency measurements of sources not identified in the NVSS, FIRST, or VLA surveys are 
given in Table 1. Unless otherwise noted, the tabulated 1.4 GHz fluxes are from NVSS. 
For the 7 sources that do not have unambiguous 1.4 GHz counterparts in NVSS, FIRST, 
or the VLA archive, we assume a 1.4 GHz flux equal to three times the survey noise. In 
reality, the true 1.4 GHz fluxes typically will be weaker and the true spectral indices of 
these sources more shallow. However, because there are only 7 sources lacking 1.4 GHz 
detections, the bias in the resulting mean spectral index for all sources in cluster fields 
computed below is small. Furthermore, because the 1.4 and 28.5 GHz measurements were 
not made contemporaneously, variability of the sources may contribute to a broadening in 
the distribution of spectral indices. 

We compute spectral indices between 1.4 and 28.5 GHz for sources in the cluster fields 
(selected at 28.5 GHz) where index a is defined by S oc v~ a . If we use only the 88 sources 
with detections in NVSS, FIRST, or the VLA maps and omit the 7 sources lacking 1.4 GHz 
detections, the resulting mean spectral index is a = 0.66 with an rms dispersion of 0.36. 
If we include limits for the 7 sources lacking 1.4 GHz detections, the mean spectral index 
a = 0.60 with an rms dispersion of 0.42, and a median index of 0.71, indicating that the 
bias due to omitting these 7 sources is small. We use the 88 sources with NVSS, FIRST, or 
VLA counterparts for the remainder of our spectral index analysis. A histogram of spectral 
indices for these 88 sources is shown in Figure 1. Characterizations of the spectral index 
distribution are given in Table 3. The distribution has a tail at low-a and is not well fit by 
a Gaussian. We therefore also compute the median of the distribution, as well as the 25th 
and 75th percentiles and find them to be 0.72, 0.51, and 0.92, respectively. 

While the beam attenuation factor is potentially a significant source of uncertainty 
for the 28.5 GHz sources, we find that excluding sources at large radii (attenuation factor 
greater than 5.0) does not change the results. We choose a maximum cutoff outer radius of 
6.6 arcmin for BIMA and 4.2 arcmin for OVRO, corresponding to a beam attenuation factor 
of about 30 and spanning a region twice the FWHM of the primary beam. When a source 
has observations from both BIMA and OVRO, we choose the one with the best combination 
of sensitivity and survey area. 

The radial distribution of spectral indices is shown in Figure 2; there is no apparent 
trend in spectral index with radius from the cluster center. We compare the spectral index 
distribution of the central regions of cluster fields (r < 0.5 arcmin) with the distribution of 
the outer regions of cluster fields (r > 0.5 arcmin) and find no significant differences. The 
mean spectral index for the central regions of cluster fields is a = 0.75 with an rms dispersion 
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of 0.24 and the mean spectral index for the outer regions of cluster fields is a = 0.63 with an 
rms dispersion of 0.38. The medians are 0.76 [0.56 (25%), 0.94 (75%)] for the central regions 
and 0.71 [0.42 (25%), 0.88 (75%)] for the outer regions. The results of a Kolmogorov-Smirnov 
(KS) test indicate that the distribution of spectral indices of sources in the inner regions of 
the cluster fields is consistent with that of the outer regions. The maximum distance between 
their cumulative distribution functions is 0.18, corresponding to an 64% probability that the 
two samples are drawn from the same distribution. Since we have few sources in non-cluster 
fields, we do not compute an average spectral index for this group. 

We compare the spectral index distribution of our m Jy cluster sources to those of some- 



what brighter field sourc es measured by other groups (IMason et al.l 120031 ; IWaldram et al. 



2003; iBolton et al.l 120041 ) ; a summary is given in Table 3. The CBI group (IMason et al. 



20031 ) finds a spectral index fr om 1.4 to 31 GHz of . 45 with an rms disp ersion of 0.37. As 
a follow-up to the 9C survey (IWaldram et al.ll2003l ). IBolton et al.l (120041 ) compute indices 
between several frequencies from 1.4 to 43 GHz, distinguishing between weak and strong 
sources. The mean spectral index for the lower fl ux sample is ~ . 4 for the l ower frequen 
ties, a nd steepens to ~ 0.9 from 15.2 to 43 GHz. IWaldram et al.l (120031 ) and IBolton et al. 



( 120041 ) also find that a greater percentage of the strong sources have a flat or rising spectrum. 



Our spectral indices are somewhat steeper than the spectral indices me asured by CBI 
and m uch steeper than those of the ~ Jy sources (a ~ 0) measured by WMAP ([Bennett et al. 



20031 ). When comparing source surveys, it is essential to consider the flux and frequency 
at which the sources are selected. We expect a survey of strong sources selected at high 
frequency to have a flatter spectral index than a survey of low flux sources selected at lower 
frequency. Our results are for a relatively low flux survey selected at high frequency, and it 
is interesting that the spectral index is relatively steep. The sources in our survey primarily 
lie in the environments of rich galaxy clusters and it is possible that there are significant 
differences between this population of sources and that found in surveys that do not target 
clusters. Our spectral indices for radio s o urces towards clusters are similar to those found 
towards clusters at 2.7 GHz by lSlee et al.l (119831 ) . They find that spectral indices are steeper 
in clusters than in the field and note a trend of shallower spectral indices with increasing 
cluster ra dius. Using our m ean spectral index of 0.66, the lower flux limit of 0.12 Jy at 2.7 
GHz from lSlee et al.l (119831 ) translates to 25 mJy at 28.5 GHz, which is slightly stronger than 
the upper limit of our source sample. 

As we look at sources with higher redshift, the emission frequency of the radiation in- 
creases. We might, therefore, expect these sources to have steeper spectral indices. However, 
the sources are selected at higher frequency which might bias the sample toward flatter spec- 
tral indices. In Figure 3, we plot spectral index as a function of cluster redshift and see no 
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clear trend. The mean and rms dispersion in the spectral index are 0.67 and 0.37 for z < 0.5, 
compared to 0.64 and 0.29 for z > 0.5, and 0.76 and 0.20 for z > 0.8. 

4. SOURCE COUNTS 
4.1. Analysis 

With the field selection effects in mind from Section 2.2, we compute the differential 
source counts as a function of flux, dN/dS, in several flux bins, accounting for the varying 
noise levels from field to field. We chose the flux bins in order to maximize the number of 
sources used and to have a similar number of sources in each bin. 

The survey boundary of each field for a given flux bin is set by the noise level of the 
field. For each flux bin and field, the minimum level in the flux bin sets the allowable beam- 
corrected noise level and the corresponding maximum attenuation radius for the field. For 
example, for > 5a sources in a flux bin of 1.5 - 2.5 mJy, the allowable beam-corrected noise 
level is 1.5/5 = 0.3 mJy. This noise level sets the attenuation radius for the field, the radius 
at which the beam attenuation factor equals the beam-corrected noise level divided by the 
uncorrected noise level. We set an outer boundary on the survey area for the field using 
the lesser of the attenuation radius or a maximum cutoff outer radius away from the field 
pointing center. We choose a maximum cutoff outer radius of 6.6 arcmin for BIMA and 4.2 
arcmin for OVRO, corresponding to a beam attenuation factor of about 30 and spanning 
a region twice the FWHM of the primary beam. We treat this as a hard maximum cutoff; 
even if the noise is sufficiently low to allow us to go to greater radii in our sampling of a 
field, we do not. The outer boundary is measured relative to the field pointing center. 

We further break the data into radial bins from the cluster center. The cluster center 
is determined by the location of the SZE decrement. For fields without a SZE decrement 
detection, the pointing center is used as the center of the field. For each field, flux bin, and 
radial bin, we compute the survey area within the boundary set by the radial bin and the 
noise level for the field. Typically the survey region for a given field, flux bin, and radial bin 
is a circle or annulus, sometimes cut off by the noise boundary. We compute the total survey 
area for each flux bin by adding up the area in all the fields. When a field has observations 
from both BIMA and OVRO, we choose the one with the best combination of sensitivity 
and survey area. 

For each field we identify all > 5a sources in the survey area that fall between the 
minimum and maximum fluxes of each flux bin. We count up the sources in each flux bin 
to get raw total source counts in the total survey area. The errors for the raw counts in 
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each bin are assumed to be Poisson distributed. Differential source counts (dN/ dS) and the 
associated errors are calculated by dividing the total raw counts in each bin by the total 
survey area for the corresponding flux bin and by the flux bin width. 



4.2. Results and Discussion 



Differential source counts (dN/ dS) , the number of sources, and the survey area for each 
flux bin are given in Figure 4 and Table 4 for the central regions of the cluster fields (r < 0.5 
arcmin), the outer regions of the cluster fields (r > 0.5 arcmin), and for the non-cluster fields. 
The error bars on dN/ dS are the Poisson errors on the raw source counts and do not include 
other sources of uncertainty. Typical raw counts of sources are ~ 4 in each flux bin for the 
inner regions of the cluster fields and ~ 8 for the outer cluster regions. We only detect two 
> 5cr sources in the 8 non-cluster fields that were selected without regard to possible radio 
source contamination. 

The differential source counts can be described by a power law, dN(S) /dS = N (S/ So)' 1 , 
where Sq — 1 mJy for this analysis. Best fits using a Markov chain algorithm that simul- 
taneously estimate the normalizations for the inner, outer, and non-cluster regions, and a 
common power-law index are shown with the data in Figure 4 and are given in Table 5. The 
best-fit common power-law index is 7 = — 1.98±0.20. Best fit power-laws for the central and 
outer cluster regions individually are also shown in Figure 4 and Table 5. As a cross-check, 
we compute dN/dS for the BIMA and OVRO fields separately and find good agreement; see 
Figure 5. All uncertainties represent 68% confidence intervals unless otherwise noted. 

Source counts are found to be greatly elevated toward the central core of the cluster 
fields. Using the normalizations from the best simultaneous fit, source counts are found to be 
a factor of 8.9l 2 *8 higher in the central regions than in the outer regions of the cluster fields. 
Counts are also elevated in the outer regions of the cluster fields relative to the non-cluster 
fields by a factor of 3.3^g. These overdensities imply that 97^5% of sources in the inner 
regions are cluster members, as are 70^37% of sources in the outer regions. A comparable 
overa bundance of radio sourc e s towa r d gala x y clusters is also s een at l ower radio frequencie s 



Rizza et all (12003 )). 



Slee et all (119831 . Il998[ ): lOwenl (119961 ): iLedlow & Owenl (119951 ): iReddv & Yunl (120041 ): 



We considered the possibility that gravitational lensing of background radio galaxies 
could produce the overabundance of detected radio sources in the direction of massive galaxy 
clusters. A gravitational lens with magnification factor /1 will modify the source counts 
to dN'(S)/dS = (dN (S / u) / dS) / fi 2 . If the unlensed source counts can be described by a 
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power law, dN(S)/dS oc S 7 , the n the sourc e counts will be changed by a factor B = 
(dN'(S)/dS)/(dN(S)/dS) = /i" 2+7 falainl 12002^ . T he mean magnification in the BIMA and 
OVRO cluster fields is estimated by ICooray et al.l (119981 ) to be /i ~ 1.4. In this analysis, we 
temporarily assume that all sources are background sources drawn from the same distribution 
and capable of being lensed. Using the best joint fit power law index of 7 = 1.98 ± 0.20 
we expect a factor of B = 0.99^^1, i.e., no change in source counts in the direction of 
the clusters due to lensing. Therefore we conclude that, regardless of the magnification, 
gravitational lensing can not be responsible for the significant excess of sources seen in the 
direction of clusters. 



In Figure 6 we compare our measurements of dN/ dS with the iDe Zotti et al.l (120051 ) 
30 GHz model and with measurements from other experiments, including WMAP, DASI, 
VSA, and CBI, which all examine non-cluster fields. We present source counts in terms of 
logio(S 5 ^ 2 dN /dS) , for ease of comparison. Counts in our non-cluster fields are consistent with 
those expected from the model and from extrapolations from other experiments, although 
with only two > 5cr sources in those fields, the sample variance is large. The source counts 
toward cluster fields have a similar power-law slope, but a higher normalization than expected 
from extrapolations of measurements of sky not concentrated on clusters. 



5. CONCLUSIONS 

From deep interferometric observations at 28.5 GHz of unresolved radio sources toward 
89 fields centered on massive galaxy clusters and 8 non-cluster fields, we find that differential 
source counts are greatly elevated in the centers of cluster fields. Counts are a factor of 
8.9+2'g higher in central regions (r < 0.5 arcmin) than in the outer regions (r > 0.5 arcmin) 
of the cluster fields. Counts in our non-cluster fields are consistent with those expected from 
models and from extrapolations from other experiments. In addition, source counts in the 
outer regions of cluster fields are a factor of 3.31^8 higher than counts in non-cluster fields. 

Using the NVSS and other surveys, we find a mean spectral index for sources in cluster 
fields between 1.4 and 28.5 GHz of a = 0.66 with rms dispersion of 0.36, where flux S oc v~ a . 
The distribution is skewed, with a median spectral index of 0.72 [0.51 (25%), 0.92 (75%)]. 
This is steeper than spectral indices of stronger (~ 20 mJy) field sources, and much steeper 
than those of much stronger (~ Jy) field sources measured by other surveys. No significant 
differences are found between the distributions of the spectral indices of sources in the inner 
and outer regions of the clusters. 

These results can be used for improving forecasts for radio source contamination of SZE 
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and CMB experiments. The cluster fields used in this work were chosen to contain massive 
clusters and the sources were identified at 28.5 GHz. We anticipate that the SZA and other 
instruments will be able to extend this work to less massive clusters and to sources identified 
at yet higher frequencies. 
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Fig. 1. — Histogram of the spectral index distribution for radio sources in cluster fields. The 
spectral index a is defined by S ~ v~ a . The distribution for 88 sources is shown with a solid 
(green) line, the distribution for the 67 sources in the outer regions of cluster fields (r > 0.5 
arcmin) is shown with a dashed (black) line, and the distribution for the 21 sources in the 
inner regions of cluster fields (r < 0.5 arcmin) is shown with a dotted (blue) line. The overall 
mean spectral index is a = 0.66 with an rms dispersion of 0.36. The mean spectral index for 
the outer regions of cluster fields is a = 0.63 with an rms dispersion of 0.38 and the mean 
spectral index for the the inner regions of cluster fields is a = 0.75 with an rms dispersion 
of 0.24. The medians are 0.72, 0.71, and 0.76 for all, outer, and inner, respectively. 
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Fig. 2. — Spectral index as a function of radius for 88 cluster sources. The figure shows no 
clear trend with radius. The half-power point for BIMA is shown as a dotted (black) line 
and the half-power point for OVRO is shown as a dashed (blue) line. Beyond the half power 
point, the beam attenuation factor becomes a potentially important source of systematic 
uncertainty. When a source has observations from both BIMA and OVRO, we choose the 
field with the best combination of sensitivity and survey area; repeat sources are not shown 
here. 
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Fig. 3. — Spectral index as a function of redshift for 88 cluster sources. The figure shows no 
clear trend with redshift. When a source has observations from both BIMA and OVRO, we 
choose the field with the best combination of sensitivity and survey area; repeat sources are 
not shown here. 
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Fig. 4. — Differential number counts, logi (dN / dS) , as a function of flux, log w (S), and best 
fit power laws for the central regions of cluster fields (r < 0.5 arcmin), outer regions of cluster 
fields (r > 0.5 arcmin), and non-cluster fields. Solid lines indicate the best fit power laws for 
each set individually and dotted lines indicate the best fits using a Markov chain algorithm 
to simultaneously estimate the normalizations and a common power law index. Using the 
best joint fit normalizations, we find that counts toward the outer regions of clusters are a 
factor of 3.3l l g higher than counts in the field. Counts toward the inner regions of clusters 
are a factor of 8.9l 2 J higher than the outer regions. The outer boundary used for the outer 
regions of cluster fields is set by the noise levels in the fields as described in Section 4.1. 
Error bars on the data come from Poisson errors on raw counts and do not include other 
sources of uncertainty. 
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Fig. 5. — Differential number counts as a function of flux for cluster fields taken at BIMA 
and OVRO. They are in good agreement. 
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Fig. 6. — Differential number counts as a function of flux. Sources are overdense in the cluster 
fields, especially within the central ar cminute . Our SZE data ar e shown along with measure - 
ments from CBI dMason et all 120031) VSA Jdearv et al.ll2005h . DASI faovac et al.ll2002h 



and W MAP Ka-band (jBennett et al.ll2003l ) as well as the 30 GHz model from lDe Zotti et al. 
(120051 ). Counts in the non-cluster fields are consistent with the model and extrapolations 
from other experiments. 



Table 1. Radio Sources in Cluster Fields 
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0.128 




















MACS J2228.5- 


f2036 


0.418 




22:28:34.4 


+20:36:37.0 


B 


0.110 


-20 


37 


1 


22:28:32.6 


+20:35:15.7 


118.6 


0.84 


0.12 


7.2 














B 


0.110 






2 


22:28:28.5 


+20:31:33.3 


346.6 


5.21 


0.72 


50.4 


MACS J2243.3 


-0935 


0.444 




22:43:21.0 


-09:35:25.0 


B 


0.618 


-35 


12 




























O 


0.121 






1 


22:43:17.8 


-09:35:08.4 


12.5 


0.93 


0.13 


j 


Aboil 2507 




0.196 


in 


22:56:51.6 


+05:30:12.2 


O 


0.072 






1 


22:56:44.0 


+ 05:31:15.5 


130.5 


8.15 


0.15 


14.4 



a No detection found in the literature. Using 3 times the NVSS noise level, we set an upper limit on the 1.4 GHz flux of 1.35 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of 0.0. 
^No detection found in the literature. Using 3 times the NVSS noise level, we set an upper limit on the 1.4 GHz flux of 1.35 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of 0.0. 
c FIRST catalog. 

d No detection found in the literature. Using 3 times the NVSS noise level, we set an upper limit on the 1.4 GHz flux of 1.35 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of 0.0. 
VLA map. Fluxes are obtained by CLEANing the images in AIPS. 

Using the 3 times the VLA map noise level, we set an upper limit on the 1.4 GHz flux of 0.45 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of -0.7. The source is detected at (TfeW by 
|Paf tridge et "al . 1 198C) with a flux of 2.59 mjy, which yields a spectral index of -0.4 from 5 to 28.5 GHz. 

s No detection found in the literature. Using the 3 times the VLA map noise level, we set an upper limit on the 1.4 GHz flux of 1.32 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of +0.2. 

h FIRST map. 

'No detection. Using the 3 times the NVSS noise level, we set an upper limit on the 1.4 GHz flux of 1.35 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of -0.4. 
J No detection. Using the 3 times the FIRST noise level, we set an upper limit on the 1.4 GHz flux of 0.45 mjy and an upper limit on the 1.4 to 28.5 GHz spectral index of -0.4. 

Note. — References: 



1. iStocke et~al~| | | 19911 ) 

2. iHughes 8z Birkinsha'wl J 19981 ) 

3. From Chandra X-ray spectrum (observations 3521 + 5010, 45 ksec exposure) we obtain z — 0.586 i 0.01. 

4. iDressler et al.| j 1 999T) 

5. iStruble 8z Roodl jl999ft 

6. iRomer et~ai"| i2000ft 

7. iGladders et all i2003l) 

8. iGioia & Luppinci I 1994) 



9. Prom Chandra X-ray spectrum (observations 3257+3582+6108, 70 ksoc exposure) we obtain z — 0.467 + 0.005. 
10. lEbeling et al.l J 199811 



12. 
13. 



LaRooue et al.l i2003h 



Edge et all <2003fl 



Crawford et al.l i 19991 ) 



Dahlc et al. 12002) 



From Chandra X-ray spectrum (observation 509, 10 kscc exposure) we obtain z — 0.452 +_ 0.005. 
Crawford et al.l i 19951 ) 



Allen et al. ( 1992) 



Luppino & Gioia 11995) 



19. From Chandra X-ray spectrum (observations 3252 + 5009, 35 ksec exposure) we obtain z — 0.48^q'qq. 

20. From Chandra X-ray spectrum (observations 3253 + 5008+5350, 35 ksec exposure) we obtain z — 0.51^^'^. 

21. iBShringer et al.l )2000l) 



23. 
24. 
25. 
26. 
27. 
28. 
29. 



Donahue et al.l i!999l ) 



Ebeling et al.l i2001bl) 



From Chandra X-ray spectrum (observations 3258 + 6110, 85 ksec exposure) we obtain z — 0.519 + 0.007. 
Gladders fc Yed j2005l) 



Schindler et al.l i l 19951 ) 



Patel et al. (2000) 



Henry et all )1997fl 



LaRoque et al.l i2006l ) 



30. From Chandra X-ray spectrum (observation 3285, 20 ksec exposure) we obtain z — 0.42^q'q2- 

31. From Chandra X-ray spectrum (observation 3260, 21 kscc exposure) we obtain z — 0.44 +_ 0.01. 



Table 2. Radio Sources in Non-Cluster Fields. 





Pointing Center 








zo.o Vjtiz oourcc 




1 A CVfr, 








rUVio 










Flux 


r Ida 


q(J2000) 


5(J2000) 


(mjyj 


O^J zUUU ) 


x( Tonnn'i 
O^JzUUU ) 


tlaams { ) r lux ^mjyj 


Unc. (mjy) 




BDF4 


00:28:04.4 


+28:23:06.0 


0.075 












BDF4 


00:28:04.4 


+28:23:06.0 


0.115 












HDF 


12:36:49.4 


+62:12:58.0 


0.084 


12:36:45.7 


+62:11:30.2 


91.5 0.48 


0.10 




BDF6 


18:21:00.0 


+59:15:00.0 


0.074 












BDF6 


18:21:00.0 


+59:15:00.0 


0.059 












BDF7 


06:58:45.0 


+55:17:00.0 


0.083 












BDF8 


00:17:30.0 


+29:00:00.0 


0.071 












BDF9 


12:50:15.0 


+56:52:30.0 


0.084 












BDF10 


18:12:37.2 


+58:32:00.0 


0.086 


18:12:16.6 


+58:29:10.1 


234.2 1.54 


0.23 




BDF11 


06:58:00.0 


+54:24:00.0 


0.090 












BDF12 


06:57:38.0 


+55:32:00.0 


0.104 












BDF13 


22:22:45.0 


+36:37:00.0 


0.097 












BDF14 


00:26:04.4 


+28:23:06.0 


0.093 












BDF15 


06:56:45.0 


+55:17:00.0 


0.076 


06:56:44.0 


+55:11:38.4 


321.7 7.07 


0.55 


4.50 


BDF16 


12:34:49.4 


+62:12:58.0 


0.092 












BDF17 


18:19:00.0 


+59:15:00.0 


0.090 












BDF18 


00:15:30.0 


+29:00:00.0 


0.087 












BDF19 


06:55:38.0 


+55:32:00.0 


0.084 












BDF20 


12:48:15.0 


+56:52:30.0 


0.088 


12:48:45.2 


+56:52:46.0 


247.7 1.36 


0.27 


4.30 


BDF21 


18:10:37.2 


+58:32:00.0 


0.083 













Table 3. Spectral indices. 



Data Set 


Frequency Range 


No. Sources 


Mean Index ± RMS 


Median Index [25%, 75%] 




Plux Limits 


Overall cluster 


1.4 to 28.5 GHz 


88 


0.66 ± 0.36 


0.72 


0.51, 0.92] 


~ 0.6 - 


10.0 mjy at 28.5 GHz (sec text) 


Inner cluster (r < 0.5 arcmin) 


1.4 to 28.5 GHz 


21 


0.75 ± 0.24 


0.76 


0.56, 0.94] 


~ 0.6 - 


10.0 mjy at 28.5 GHz (sec text) 


Outer cluster (r > 0.5 arcmin) 


1.4 to 28.5 GHz 


67 


0.63 ± 0.38 


0.71 


0.42, 0.88] 


~ 0.7 - 


8.0 mjy at 28.5 GHz (see text) 


Mason et al. 2003 (CBI) 


1.4 to 31 GHz 


56 


0.45 ± 0.37 








21 mjy at 31 GHz 


Bolton ct al. 2004 (9C follow-up) 


1.4 to 4.8 GHz 


124 




0.44 


0.05, 0.76] 




25 mjy at 15.2 GHz 


Bolton ct al. 2004 (9C follow-up) 


4.8 to 15.2 GHz 


124 




0.39 


0.06, 0.95] 




25 mjy at 15.2 GHz 


Bolton ct al. 2004 (9C follow-up) 


15.2 to 43 GHz 


124 




0.87 


0.42, 1.20] 




25 mjy at 15.2 GHz 


Bolton ct al. 2004 (9C follow-up) 


1.4 to 4.8 GHz 


70 




0.24 [ 


-0.12, 0.64] 




60 mjy at 15.2 GHz 


Bolton ct al. 2004 (9C follow-up) 


4.8 to 15.2 GHz 


70 




0.27 


0.02, 0.70] 




60 mjy at 15.2 GHz 


Bolton ct al. 2004 (9C follow-up) 


15.2 to 43 GHz 


70 




0.67 


0.38, 1.03] 




60 mjy at 15.2 GHz 



I 

CO 

o 

I 



Table 4. Number Counts for > 5a Sources at 28.5 GHz 



loglO(S) (mjy) 


radius range (arcmin) a 


raw counts 


area (arcmin 2 ) 


logio(dN/dS) (arcmin 2 mjy 1 ) 


field type 


-0.15 to 0.00 
0.00 : 0.15 
0.15 : 0.30 
0.30 : 0.45 
0.45 : 0.60 
0.60 : 0.75 
0.75 : 0.90 


> 0.5 

> 0.5 

> 0.5 

> 0.5 

> 0.5 

> 0.5 

> 0.5 


10 

6 

10 
6 
11 

8 
7 


783 
1541 
2420 
3306 
4217 
5104 
5951 


1 oc: + 0.15 

± - JO -0.l6 
9 n o+0.20 
z.uo_ 22 

r, i t-+0.15 

- 2 - 15 -0.15 
-2 66+ ' 20 

Z.DD_g 22 

-2 65 +0 15 

-3 02+° 17 
_o 2 „+0.l9 


cluster 
cluster 
cluster 
cluster 
cluster 
cluster 
cluster 


-0.20 to 0.10 
0.10 : 0.40 
0.40 : 0.70 
0.70 : 1.00 


< 0.5 

< 0.5 

< 0.5 

< 0.5 


4 
6 
3 
2 


42.0 
63.4 
69.0 
69.5 


-0.82t££ 

-1 19 +0.20 
± - lz -0.22 
-1 76+ ' 30 

-2 24+ ' 37 


cluster 
cluster 
cluster 
cluster 


0.13 : 
0.50 : 


0.50 
0.87 


> 0.0 

> 0.0 


1 
1 


438 
724 


o 4 q+0.52 


non-cluster 
non-cluster 



a Max radius varies for individual fields since set by noise. See text. 



Table 5. Power-law Fits to dN/dS 



Field type 


Best indiv. fit index 


Best indiv. 


fit normalization 


Best joint fit index 


Best joint fit normalization 


Inner cluster (r < 0.5 arcmin) 
Outer cluster (r > 0.5 arcmin) 
Non-cluster 


-1.67 ± 0.37 
-2.02 ± 0.22 


132l 4 J 8 x 10" 
20.9±4 3 x 10 


~ 3 mjy - 1 arcmin -2 
~ 3 mjy -1 arcmin -2 


-1.98 ±0.20 
-1.98 ±0.20 
-1.98 ±0.20 


174^45 X 10~ 3 mjy" 1 arcmin^ 2 
19.5lg;° X 10~ 3 mjy" 1 arcmin" 2 
5.90±|;™ X 10~ 3 mjy" 1 arcmin- 2 



